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ABSTRACT 

The  problem  of  determining  whether  a  line  of  bearing 
measured  by  a  local  surface-based  sensor  coincides  with  a 
threat  whose  position  has  been  previously  estimated  by  an 
ocean  surveillance  sensor  is  addressed.   Uncertainties  in 
the  position  estimate  of  the  threat,  in  the  bearing  error 
and  in  the  position  estimate  of  the  sensor  are  considered 
in  measuring  the  probability  that  the  threat  lies  on  a  given 
bearing  from  the  sensor.   A  TI-59  calculator  program  is 
developed  which  calculates  this  likelihood  when  the  threat 
location  density  can  be  assumed  to  be  bivariate  normal. 
Computations  required  when  significant  time  has  elapsed 
since  the  original  estimate  of  threat  location  when  the  den- 
sity can  no  longer  be  considered  bivariate  normal  are 
discussed. 
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I.   INTRODUCTION 

As  long  range  surface  launched  weapons  systems  continue 
to  be  introduced  into  the  fleet,  the  operational  commander 
is  increasingly  faced  with  the  problem  of  being  able  to 
launch  weapons  at  targets  located  beyond  the  horizon. 
Successful  employment  of  such  weapons  is  dependent  not  only 
on  the  ability  to  detect,  classify  and  localize  targets  at 
considerable  distances,  but  also  on  the  ability  to  distinguish 
the  true  target  from  a  potentially  larger  field  of  false 
targets.   While  long  range  ocean  surveillance  sensors  may  be 
of  assistance  in  the  identification  and  localization  of   tar- 
gets, the  information  provided  may  not  be  refined  sufficiently 
to  permit  effective  targeting  of  long  range  weapons  on  that 
basis  alone.   The  on-scene  commander  must  in  general  rely  on 
additional  data  on  target  location  gathered  locally  and  close 
to  the  time  of  weapons  launch  for  accurate  targeting.   Thus, 
he  must  still  be  able  to  detect  and  track  the  desired  target 
and  be  able  to  distinguish  it  from  other  targets  within 
range  of  his  sensors. 

The  procedures  developed  in  this  paper  are  designed  to  be 
of  assistance  in  addressing  the  last  of  these  problems.   They 
are  applicable  when  the  information  available  is  an  error 
ellipse  around  a  threat  location  estimated  by  an  ocean  sur- 
veillance sensor  and  bearings  only  data  generated  by  a  local 
surface-based  sensor.   The  question  of  whether  a  target 


detected  by  the  local  sensor  is  the  same  as  that  whose 
estimated  position  was  provided  by  an  external  sensor  can 
only  be  addressed  if  information  is  available  on  the  loca- 
tions and  tracks  of  all  possible  targets  within  range  of  the 
local  sensor.   Since  such  data  is  generally  not  available, 
this  paper  does  not  attempt  to  answer  that  question,  but 
rather  develops  a  method  by  which  bearing  information  from 
different  sensors  can  be  compared  as  to  the  likelihood  of 
each  bearing  being  associated  with  the  threat  identified 
previously.   It  is  envisioned  that  these  likelihoods  can 
be  then  used  to  induce  an  ordering  among  bearing  data  gathered 
by  different  sensors  or,  conceivably,  conflicting  data 
gathered  by  one  sensor.   The  ordering  would  be  based  on  the 
likelihood  that  each  bearing  will  contribute  to  refining  the 
original  estimate  of  the  location  of  the  target  of  interest. 
This  information  could  be  applied  in  a  number  of  ways:   as 
a  guide  to  allocation  of  more  capable  sensor  resources  for 
purposes  of  obtaining  targeting  information;  as  a  guide  for 
allocation  of  weapons  against  more  than  one  threat;  as  a 
means  of  pre-processing  data  before  entering  it  into  a  target 
motion  model,  thereby  reducing  the  chance  of  introducing 
unrelated  data. 

To  determine  the  likelihood  that  a  given  line  of  bearing 
and  the  threat  coincide,  consideration  was  given  to  the 
uncertainties  inherent  in  estimation  of  target  position, 
in  the  measurement  of  bearings  by  a  particular  sensor  and 
in  estimation  of  sensor  location.   It  is  assumed  that  at  some 


time  t  ,    an  ocean  surveillance  sensor  detects  a  threat  whose 
o 

position  is  estimated  to  be  within  an  elliptical  region  with 
p,  X  100%  certainty.   The  estimated  position  data  are  received 
and  converted  by  the  on-scene  commander  into  a  probability 
distribution  described  by  a  truncated  bivariate  normal  density 
function. 

It  is  further  assumed  that  the  standard  error  a.  charac- 

3 

teristic  of  the  local  sensor  is  known.   The  sensor  bearing  3 

with  bearing  error  c^n  is  then  projected  from  the  sensor 

position  through  the  threat  density  function. 

Since  sensor  position  relative  to  the  target  may  itself 

be  subject  to  navigation  error,  the  uncertainty  is  introduced 

as  a  truncated  circular  bivariate  normal  distribution  centered 

at  location  (u  ,v    )    with  standard  deviation  a. 
o   o 

A  TI-59  calculator  program  is  developed  which  estimates 
the  likelihood  that  the  threat  identified  by  an  external 
sensor  lies  along  bearing  S,  given  the  threat  distribution, 
the  bearing  error,  and  sensor  position  distribution  relative 
to  the  threat. 

The  theoretical  basis  for  this  calculation  is  presented 
in  Chapter  II.   The  algorithms  used  in  designing  the  calcu- 
lator program  are  described  in  Chapter  III.   A  program  listing 
and  verbal  flow  are  provided  in  Appendix  A  along  with  instruc- 
tions for  the  user.   Appendix  B  contains  a  verbal  flow  of  a 
program  designed  for  use  when  considerable  time  has  elapsed 
since  the  initial  estimate  of  the  location  of  a  moving 
target. 


II.   THEORETICAL  BASES 

The  general  approach  to  determining  the  likelihood  that 
a  measured  line  of  bearing  B  is  the  true  bearing  from  the 
sensor  to  the  threat  identified  and  localized  by  an  external 
sensor  is  discussed  in  this  chapter.   Calculations  required 
when  using  threat  position  information  both  as  initially 
generated  by  the  external  sensor  at  time  zero,  t  ,  and  as 
distorted  to  account  for  an  intervening  time  late,  t  ,  are 
discussed.   Uncertainties  in  bearing  measurement  and  sensor 
position  are  included. 

Initially,  assume  that  sensor  position  is  known  with 

certainty.   Let  3_  be  the  true  bearing  of  the  threat  from  the 

sensor.   Since  threat  location  is  uncertain,  6„  is  a  random 

quantity  with  probability  density  function  f      (6m)-   Let 

^TRUE 
f„(S;6m)d6  be  the  probability  that  the  errors  in  bearing 

measurement  are  such  as  to  give  rise  to  a  bearing  on  the 

threat  in  the  interval  d6  about  the  observed  value  s  when 

the  true  bearing  is  3  .   Then  the  likelihood  of  observing  a 

bearing  3  is: 


)   =     /    f«(S;6m)  f«     (3m)  dBm  .    (1) 
all  6^   ^     ^    ^TRUE   ^     ^ 


The  density  f      (3  )  is  determined  by  the  probability 

^TRUE 
density  function  of  threat  location  at  time  t,  f^   (x,y;t; 

A,  1 
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Given  the  probability  density  function  of  threat  location 
at  time  t  and  assuming  that  bearing  errors  are  normally  dis- 
tributed with  mean  zero,  consideration  is  limited  to  the 
probability  that  the  observed  bearing  is  the  true  bearing 
of  the  threat,  given  the  threat  is  contained  in  a  planar 
region  E  and  the  true  bearing  B   is  in  an  interval  about  g 
with  upper  and  lower  bounds  B.t  and  B    .      The  region  E  is 

U  Li 

selected  to  be  the  minimum  area  planar  region  which  contains 
the  threat  with  a  specified  high  probability  p, .   The  inter- 
val (0  ,8^)  is  selected  so  that,  for  all  the  possible  values 
of  6_  contained  in  the  interval,  B  is  within  an  interval  of 
specified  high  probability  p^  around  Bm-   The  fan  (3,8.,) 

Z  1  Li    U 

is  symmetric  about  B  so  that  it  represents  the  minimum  area 
region  which  meets  the  above  criterion.   The  likelihood  of 
observing  B  when  the  threat  is  in  E  and  Bm  is  in  (6-,./ Br,)  is 
determined  by  integrating  expression  (1)  above,  over  a  region 
defined  in  the  manner  of  the  shaded  area  of  Figure  1(a). 
If  sensor  position  is  not  certain,  assume  that  it  is 
distributed  in  accordance  with  a  circular  bivariate  normal 
distribution  f    (u,v) .   Again,  consideration  is  limited  to 
determining  the  likelihood  of  observing  a  bearing  B  given 
that  the  threat  is  contained  in  a  planar  region,  E,  6m  is  in 
(6t/3tj)/  and  the  sensor  is  contained  in  a  planar  region  C. 
The  region  C  is  selected  as  the  minimum  area  planar  region 
which  contains  the  sensor  with  a  high  probability  p...   The 
region  of  E  over  which  f ' (B)  is  evaluated  expands  as  the 
shaded  region  of  Figure  1(b). 
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%'V 


(a)   Sensor  Position  Certain 


(b)   Sensor  Position  Uncertain 


FIGURE  1.   REGION  OF  INTEGRATION 
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since  the  measurement  errors  involved  in  estimating  the 
threat  position,  the  sensor  position  and  the  bearing  angle 
arise  from  different  measurement  procedures,  it  is  reason- 
able to  assume  that  errors  are  independent.   The  probability 
densities  required  to  perform  the  above  calculations  can  be 
estimated  as  follows. 

At  time  t   an  ocean  surveillance  sensor  estimates  the 
o 

position  of  the  threat  to  be  located  within  an  elliptical 
area  characterized  by  the  parameter  set  E  =  {X,Y,e,A,B}  with 
confidence  p,  x  100%.   The  elements  of  the  parameter  set  E 
are  respectively:   X  the  latitude  of  the  estimated  threat 
position,  Y  the  longitude  of  estimated  threat  position, 
6  the  orientation  of  the  major  axis  from  true  North,  A  the 
length  of  the  semi-major  axis,  and  B  the  length  of  the  semi- 
minor  axis.   Since  the  measurement  errors  in  determining 
threat  position  are  generally  assumed  to  be  normally  distributed, 
the  ellipse  characterized  by  E  represents  the  minimal  area 
p,  X  100%  confidence  region  about  the  mean  (X,Y) .   Treating 
this  ellipse  as  the  p,  probability  region  of  a  bivariate 
normal  distribution,  a  density  function  for  the  threat  posi- 
tion at  time  t   can  be  estimated.   For  convenience,  locate 
o 

the  origin  of  a  rectangular  coordinate  system  at  the  center 

of  the  ellipse,  (X,Y) ,  with  positive  x-axis  located  along 

the  major-axis  of  the  ellipse  at  a  bearing  9  from  true  North. 

Assume  a  flat  earth  in  the  region  of  interest.   Let  t  =0. 

^  o 

The  mean  of  the  threat  position  density  f    (x,y;0)  is  then 

A./  X 

the  point  (0,0).   The  variances  in  the  X  and  Y  directions  can 
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be  derived  from  the  fact  that  the  region  with  minimal  area 
which  contains  the  threat  with  probability  p,  is  a  k-sigma 
ellipse  where  k  is  determined  from  the  relationship: 


-k2/2 
P  [threat  located  in  k-sigma  ellipse]   =   1  -  e   ^ 


Thus, 


a^      =   (A/k)  ^ 


and 


a^^   =   (B/k)^ 


where 


k   =   /-2  ln(l  -  p^) 


X  and  Y  are  assijmed  to  be  independent. 

If  the  course  and  speed  of  the  threat  are  known  with 
certainty  to  be  ^   and  s  respectively,  the  probability  density 
of  the  threat  position  at  time  late  t-  can  be  shown  to  be 
again  a  bivariate  normal  with  mean  (st  cos  ( 6-4^)  ,  st  sin  (8-4^)  ) 

Li  Li 

2     2 

and  variances  o      ,    a      ,      The  p,  probability  region  of  the 

X      Y  L 

density  at  time  t  would  then  be  an  ellipse  congruent  to  that 

Li 

characterized  by  the  set  E  above  but  centered  at  the  point 
(st^  cos  (9-iM /St_  sin(e-i|;)  )  (Figure  2)  [Ref.  1]. 

Li  Li 
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T 


;t^sin(9-4^)) 


FIGURE  2.   TIME  LATE  ELLIPSE,  KNOWN  COURSE  AND  SPEED 
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In  some  cases,  the  motion  of  a  submarine  on  patrol  in  a 
large  area  can  be  characterized,  when  t_  is  large,  by  an 
expansion  of  the  probability  area  with  time  at  a  rate  D. 
The  result  in  such  a  case  is  that  f„  ^{yL,Y;t^)    is  still 

A  ,  X  Li 

2 
bivariate  normal  with  mean  (0,0)  but  with  variances  a^  +  Dt^. 

X  Li 

2 

and  a   +  Dt  . 

Y  Li 

When  the  motion  of  the  threat  cannot  be  described  in  the 
above  manner,  and  the  course  and  speed  are  not  known,  but 
assumed  to  be  distributed  according  to  the  densities  f-^iii) 
and  fc,(s),  determining  the  probability  density  f„  ^(x,y;tT.) 

O  A  ,  X  Li 

is  a  considerably  more  complex  problem.   Let  (x  (s,ij;),y  {St'\i)) 
be  the  coordinates  of  the  point  at  which  the  threat  would 
have  to  be  located  at  time  zero  in  order  to  reach  the  point 
(x  , y  )  at  time  t^.  if  the  threat  speed  were  s  and  course  ^. 

Li    Li  Li 

Thus,  X  {s,i))    =   X  -St..  cos  (6-4;)  and  y^(s,'4;)  =  y  -st  sinCe-ij;) 

O  Li      Li  O  Li      Li 

(Figure  3) .   Then, 


360 


f^(.i;)fg(s)  d^ds    [Ref.  1] 


This  density  is  no  longer  normal.   In  the  special  case  when 

2      2 
A  =  B,  i.e.,  a„  =   o^    ,    ^   has  a  uniform  distribution  over 

X       X 

the  interval  (0*,360°)  and  S  is  known  with  certainty,  as 
shown  in  Reference  2,  the  density  changes  with  time  as  in 
Figiire  4. 
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\'h 


O   O 


FIGURE  3.   TIME  LATE  POSITION,  COURSE  AND  SPEED  UNCERTAIN 
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t=0 


t=2 


t=3 


FIGURE  4.   TIME  LATE  THREAT  DENSITY,  SPEED  KNOWN 
COURSE  UNIFORM  (O^^SSO") 
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The  distribution  of  the  bearing  error  measured  by  the 

local  sensor  can  be  estimated  if  the  standard  error  of  the 

sensor,  a    ,    is  known.   The  bearing  error  is  then  assumed  to 
p 

2 
be  normally  distributed  with  mean  zero  and  variance  a^    . 

p 

Assume  that  the  estimated  sensor  location  (u  ,v  )  is 

o   o 

accurate  to  within  R  nautical  miles  with  p-  x  100%  confidence. 
The  density  f    (u,v)  of  sensor  location  can  be  assumed  to  be 

U  /  V 

a   circular   bivariate  normal  with  mean  (u  ,v  )  and  variance 

o   o 

2     2   2 

a   =  R  /k  .   The  value  of  k  is  determined  from  the  relation- 
ship: 


2 

-k  /2 
P[(u,v)  contained  in  k-sigma  region]   =   1  -  e 


where  the  probability  on  the  left  is  p_  in  this  case. 

Since  the  evaluation  of  f'(6)  considering  sensor  position 
density  and  bearing  error  density  does  not  generalize  to  a 
closed  form,  algorithms  are  developed  in  the  remainder  of  this 
paper  for  estimation  of  the  probability  (likelihood)  that  the  threat 
identified  by  an  external  sensor  lies  on  a  line  of  bearing 
measured  by  a  local  sensor  given  that  the  threat  is  in  region 
E,  the  sensor  is  in  region  C  and  the  true  line  of  bearing 
lies  within  {&,&). 
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III.   ALGORITHMS 

As  indicated  in  the  previous  chapter,  the  variety  of 
geometrical  situations  which  can  arise  depending  on  the 
location  of  the  sensor  relative  to  the  estimated  threat  loca- 
tion precluded  development  of  a  generalized  analytical  pro- 
cedure.  Rather,  algorithms  are  developed  in  this  paper  for 
numerically  evaluating  f'(3).   The  procedures  applicable  to 
a  bivariate  normal  threat  location  density  have  been  imple- 
mented on  a  TI-59  calculator.   Appendix  A  contains  a  listing 
of  that  program.   The  calculations  required  when  the  time 
late  threat  location  density  is  no  longer  normal  exceeded 
the  available  program  capacity  of  the  TI-59  and  therefore 
have  not  been  implemented.   A  detailed  verbal  flow  is  provided 
at  Appendix  B  for  future  implementation  on  a  larger  machine. 
The  algorithms  used  in  both  situations  are  described  in  this 
chapter. 

A.   BIVARIATE  NORMAL  THREAT  LOCATION  DENSITY 

This  case  includes  situations  (1)  where  the  time  elapsed 
since  generation  of  the  initial  error  ellipse  by  the  ocean 
surveillance  sensor  is  negligible,  (2)  where  the  motion  of 
the  threat  can  be  assumed  to  be  random  in  the  manner  described 
above,  and  (3)  where  course  and  speed  of  the  threat  are 
assumed  to  be  known  with  certainty.   With  appropriate  modi- 
fications to  the  input  data,  all  three  of  these  situations 
can  be  addressed  using  the  program  contained  in  Appendix  A. 
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In  situation  (1)  the  data  entered  are  the  parameters  of  the 
ellipse  as  generated  by  the  ocean  surveillance  sensor.   In 
(2) ,  the  location  and  orientation  from  North  of  the  ellipse 
is  the  same  as  originally  generated,  but  the  size  of  the 
ellipse  expands  at  some  constant  rate  of  area  per  unit  time 
which  must  be  estimated  by  the  user.   This  rate  D,  times  the 
elapsed  time,  t  ,  yields  the  factor  which  must  be  added  to 

Li 

the  semi-^najor  and  semi-minor  axes  of  the  original  ellipse. 
That  is,  if  the  original  error  ellipse  is  a  p  ^  100%  confi- 
dence ellipse,  the  semi-major  and  semi-minor  axes  of  the 
diffused  ellipse  will  be  input  as  A'  and  B'  respectively: 


A'   =  Va^  +  (-2  ln(l-pj_))Dt^ 

B'   =  VB^  +  (-2  ln(l-p^))Dt 

In  situation  (3)  the  dimensions  and  orientation  of  the  time 
late  ellipse  are  the  same  as  those  of  the  original  error 
ellipse,  but  the  center  of  the  ellipse  is  displaced  from 
its  original  position  by  the  known  velocity  times  elapsed 
time.   The  updated  position  of  the  error  ellipse  is  treated 
as  the  origin  of  the  rectangular  coordinate  system  for  this 
situation  and  all  linear  measurements  are  made  relative  to 
this  system.   All  angular  measurements  are  made  from  true 
North. 

Estimation  of  the  likelihood  that  the  threat  lies  on  a 
bearing  from  the  local  sensor  given  that  the  threat  is  located 
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in  a  p,  X  100%  confidence  ellipse  and  the  true  bearing  lies 
within  the  bounds  {&    ,Q    )    with  confidence  p^  x  100%  proceeds 

Li    U  Z 

as  follows; 

1.  Estimate  parameters  of  the  bivariate  normal  density 

2     2 

f^  Y^^'^^  °^  threat  location:  u^   =    My   =    0,  ^^^y  /  ^Y  * 

2.  Determine  3-  and  Q      such  that  an  interval  of  length 

Li  U 

2k. o.    centered  on  either  6_  or  3,7  would  contain  B,  the 

p  p  L       U 

measured  bearing,  with  probability  p_:   3,.  =  3  -  k  a,, 

^        Li  3   3 

3.  Determine  sensor  location  coordinates  relative  to 
the  origin  of  the  threat  ellipse. 

4.  Subdivide  the  angular  interval  (3,3)  into  2n  sub- 
intervals.   Each  subinterval  k  intersects  the  ellipse  in  a 
strip  with  average  width  W,  which  corresponds  to  A3x.« 

5.  At  the  midpoint  of  each  subinterval  k  determine  the 
equation  of  the  line  through  the  sensor  position  at  the  true 

bearing  3,  from  North. 

Y-c 

6.  Let  the  equation  of  the  line  of  bearing  6,  be  X  =  . 

km 

Then  the  plane  perpendicular  to  the  xy-plane  which  contains 
this  line  intersects  the  bivariate  normal  threat  density  in 
a  curve  whose  equation  is 

2^  2   2       2' 

g(y)   =  T-^ e      "^      ^  (2) 

Y-c 
found  by  making  the  substitution  X  =  — —  m  the  density 

function  f^   (x,y) .   It  will  prove  convenient  to  expand  the 

right  side  of  (2)  as  follows  [Ref.  3]: 
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_1  c 

"2  2^22 

g(y)   =  "^                 e  ^      ^          (3) 

^               ^  2.22            2 


2^2   2       2^2   2   2^^^"   2^2   2^ 

[- e  ] 

/2T  a  am 


,/^ 


For  computational  purposes,  assume  that  the  width  W, 
of  the  region  of  the  ellipse  cut  out  by  the  angular  subinter- 
val  around  B,  will  be  nearly  constant  through  the  ellipse. 
Let  the  points  of  intersection  of  the  line  of  bearing  6, 
with  the  threat  ellipse  be  (^1,1'^,,)  and  (^V2'^k2^*   Then 
approximate  the  volume  of  the  normal  density  over  this  region 
by  the  absolute  value  of  the  product  of  the  area  under  the 
curve  g  between  Y,  ,  and  Y,  _  and  W,  (Figure  5)  .  Observe  that 
the  term  in  brackets  in  equation  (3)  above  is  the  density 
function  of  a  univariate  normal  random  variable  with  mean 

2 


and  variance 


2^22 

o^  X 


2   2   2 

m  a^  a^ 

T~     2  2  • 

o^  X 


The  term  preceding  the  brackets  in  (3)  is  the  slope  of  the 

line  of  bearing,  m,  times  the  density  of  a  univariate  normal 

2    2   2 

random  variable  with  mean  zero  and  variance  Oy    +  m  a^ 
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FIGURE    5 
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evaluated  at  c.   Thus,  evaluating  the  expression 

J 
Y, 


k2 

\\  ^^    S        g(y)  dy 

"kl 


is  equivalent  to  the  computation 


Wj^m(      -"      4)(  -))  (^(Z^)  -  $(Z^)) 


y   2~~2   2  J     T~~2       2 


where 


2       /  2  ^  2   2 

ca^       i/a^  +ni  <^x 

^i   "    ^^ki        2~T2~~2'^  ^     E^T^;;^        ^ 

a   +  m  a^      ,  X  Y 


for  i  =  1,2,   ^  is  the   P[Z_<_z]  when  Z  is  a  standard  normal 
random  variable,  and  ^)  is  the  density  function  of  a  standard 
normal  random  variable. 

If  the  substitution  Y  =  mX  +  c  were  made  in  the  density 
f    (x,y),  an  analogous  situation  would  arise  with  the  limits 
of  the  integration  being  X,  ,  and  X,  „. 

7.   If  the  value  of  the  density  under  the  bivariate 
normal  curve  over  the  region  of  the  ellipse  defined  by  the 
subinterval  A6^  is  weighted  by  the  probability  of  observing 
a  bearing  error  (6,-6)  the  result  is  the  probability  of 
observing  the  bearing  of  the  threat  as  6  when  the  true  threat 
location  is  in  the  segment  of  the  ellipse  defined  by  A6i^. 
This  probability  is  approximated  as  follows: 
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^2 


Since  the  bearing  error  is  assumed  normally  distributed  with 

2 

mean  zero  and  variance  a      ,    the  value  of  f  (3;6, )  can  be 

p  B      K 


determined  by  the  expression 


where  cj)  is  as  above  the  density  function  of  a  standard 
normal  random  variable. 


f(3)   =     I    fo(e;B.)  f,     (3.)  d8, 
all  6^   ^     ^    ^TRUE   ^     ^ 


is  approximated  by 


Y 
2n  k2 

^"-^   ^ki 


The  value  of  the  sum  is  determined  by  repeating  steps 
A. 6  and  A. 7  above  at  the  midpoints  of  each  of  the  2n  sub- 
intervals  defined  in  step  A. 4  and  summing  the  results  of 
each  of  these  calculations.   Obviously,  the  finer  the  sub- 
division of  (6T/3rf),  the  more  accurate  will  be  the  estimate 
J-i   u 

of  the  likelihood,  but  also  the  longer  the  calculation  will 
take. 
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B.   SENSOR  POSITION  UNCERTAINTY 

The  result  of  the  above  calculation  will  be  the  likeli- 
hood f (3)  that  the  threat  lies  along  the  bearing  measured 
by  the  sensor  given  that  the  threat  is  located  within  the 
threat  ellipse  and  the  true  bearing  is  within  the  interval 
(6  , S„)  and  the  sensor  is  at  the  position  used  to  perform 
the  calculation.   We  next  will  introduce  additional  calcu- 
lations that  are  required  to  account  for  the  fact  that  the 
sensor  position  is  not  known  with  certainty. 

1.  If  the  position  of  the  sensor  is  estimated  as  being 
within  R  nautical  miles  of  (^-./V  ),  its  assumed  coordinates 
in  the  xy-system,  with  p^  x  100%  confidence,  estimate  the 

parameters  of  the  sensor  location  density  f    (u,v)  with 

2 

mean  zero  and  variance  a  :   mean  =  (u  /V  ) ,  variance, 

2     2 
a   =  R  /(-2  In(l-p-))  in  threat  centered  coordinates. 

2.  The  bearing  measured  by  the  sensor  is  B  regardless 

of  the  sensor  location.   Assume  that  the  area  of  intersection 
of  the  angular  wedge  (6y/3u)  ^^^  the  threat  ellipse  does  not 
change  significantly  as  the  sensor  position  is  moved  along 
the  line  of  bearing  g.   Then  the  effect  of  the  bivariate 
normal  distribution  of  sensor  location  can  be  approximated 
by  considering  only  the  univariate  normal  density  along  a 
line  through  (u  ,v  )  perpendicular  to  3.   Repeat  the  calcu- 
lations in  steps  A. 3  through  A. 8  above  with  the  sensor  located 

at  each  of  the  three  points  (u  ,v  ),  (u  +. 97acos ( e-B-90) , 

o   o     o 

V  +.97asin(e-B-90) ) ,  and  (u  -. 97acos (6-3-90) ,v  -97asin ( 9-3-90) ) 
o  o  o 

If  the  line  through  (u  ,v    )    perpendicular  to  3  is  subdivided 


27 


symmetrically  about  (u  ,v  )  such  that  1/3  of  the  univariate 

o   o 

normal  density  lies  above  each  subinterval,  the  three  points 
chosen  above  represent  the  "center  of  gravity"  of  each  third 
of  the  density  (Figure  6) . 

3.  If  R  is  chosen  to  include  a  significant  proportion 
of  the  sensor  density,  i.e.,  on  the  order  of  2a  or  greater, 
the  probability  of  the  sensor  being  located  in  each  of  the 
three  regions  is  approximately  1/3.   Thus,  if  p^  is  on  the 
order  of  .86,  multiply  each  result  in  step  B.2  by  1/3. 

4.  Summing  the  results  of  steps  B.2  and  B.3  yields  an 
estimate  of  the  likelihood  that  the  threat  lies  at  bearing 
6  given  the  threat  is  in  the  p,  x  100%  confidence  ellipse, 
6_  is  in  is,    ,  Q    )    and  the  sensor  is  in  a  p^  ^  100%  confidence 
region.   That  is,  the  likelihood  is  estimated  by 


^   1 

f '  (3)   =    [   f  f .  •  (B)  , 


3  =  1  '    ^ 


where  j  is  the  index  of  sensor  position  in  figures  5  and  6. 

Instructions  for  application  of  the  TI-59  program  to 
calculate  the  above  are  included  in  Appendix  A. 

C.   THREAT  DISTRIBUTION  NOT  BIVARIATE  NORMAL  AFTER  TIME 
LATE  ELAPSED 

The  basic  approach  to  evaluating  f'(3)  when  the  time  late 

distribution  of  the  threat  is  no  longer  bivariate  normal 

is  the  same  as  that  just  discussed.   The  principal  difference 
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I .86a-4 


(u  -.97aoos(6-B-90) , 

V  -.97asin(9-6-90)) 
o 

=  (U3,V3) 


(u  +.97ocos(e-S-90) , 
o 

V  +.97  sin{9-6-90)) 
o 


=  (u^^v^) 


FIGURE  6.   ESTIMATE  OF  SENSOR  POSITION  DENSITY 
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arises  from  the  fact  that  the  time  late  density  is  signi- 
ficantly more  complex  in  this  case. 

The  threat  distribution  becomes  distorted  from  the  normal 
after  some  time  late  when  the  course  and  speed  of  the  threat 
are  constant  but  not  known  with  certainty.   Application  of  the 
method  described  herein  requires  that  the  user  assume  a  dis- 
crete distribution  of  the  speed  of  the  threat  with  upper  and 
lower  bounds  Sj.  and  S, ,  respectively.   In  addition,  the  threat 
course  is  assumed  to  be  uniformly  between  0°    and  360°. 

The  density  of  threat  location  after  some  time  late  t 
when  the  threat  speed  is  s.  then  becomes  [Ref.  1]: 


2 

,      360     ,  x-s  .  t  cos  (  9-iij)  ) 

f^   (x,y;t^,s.)   =   27^;^/    eKp[-^{ i ^ 

X  Y  0  a 

A 


(y-s.t^  sin(e-'J;)  )     ,, 
+  ^^ )  ]  ^^ 


2  '  ■■  360 

"y 

where  6  is  the  bearing  of  the  major  axis  from  North.  Note 
that  the  new  threat  density  is  still  centered  at  the  same 
position  as  the  time  zero  ellipse  but  its  shape  changes  as 
in  Figure  4  of  Chapter  II.  If  the  threat  speed  is  s.,  and 
the  course  is  uniformly  distributed  over  (0'',360'')/  the 
outer  limit  of  the  new  planar  region  containing  the  threat 
after  time  t   has  elapsed,  given  that  it  was  originally 

Li 

located  in  the  p,  x  100%  confidence  ellipse  with  semi-axes  A 
and  B,  can  be  represented  by  an  ellipse  with  semi-major  axis 
A+s .  t^  and  semi-minor  axis  B+s-t^.  (Figure  7).   Thus  the 

1   Li  1   Li 
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FIGURE  7.   TIME  LATE  PLANAR  REGION 
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region  over  which  the  density  will  be  evaluated  is  still 
the  intersection  of  an  ellipse  with  an  angular  wedge, 
recalling  however  that  the  density  is  no  longer  normal. 
Further,  the  new  elliptical  region  does  not  represent  a 
p,  X  100%  confidence  region  of  the  time  late  density. 

Calculation  of  the  likelihood  that  the  threat  lies  along 
bearing  6  given  that  the  threat  lies  in  the  p,  x  100%  ellipse 
at  time  zero,  that  the  true  bearing  lies  in  {Q    ,Q.),    that 

J-i   u 

the  sensor  lies  in  the  p-  x  100%  circle  and  that  the  speed  is 
s.  proceeds  as  follows: 

1.  Estimate  the  parameters  of  the  original  normal  dis- 

2     2 
tribution:   mean  =  (0,0),  a^  =  A  /(-2  ln(l-p,)), 

a^^  =  B^/(-2  ln(l  -  p^) )  . 

2 .  Determine  the  upper  and  lower  bounds  on  the  true 

bearing  wedge,  6  +  k^a^  and  3  -  k.a^ . 

p  p         p  p 

3.  Determine  the  time  late  planar  region  as  the  ellipse 
with  semi-major  axis  equal  to  A + s . t^  and  semi -minor  axis 
equal  to  B  +  s . t  . 

4.  Determine  the  position  of  the  sensor  relative  to  the 
ellipse  center. 

5.  Subdivide  the  bearing  fan  {Q  -  k.a ^ , B  +  k-C .)    into 

p  p      p  p 

2n  subintervals. 

6.  At  the  midpoint  of  each  subinterval,  determine  the 
equation  of  the  line  through  the  ship  position  at  that 
bearing  B  . 

7.  Let  the  equation  of  the  line  of  bearing  6,  be 

Y-c 

X  =  .   Then,  the  plane  perpendicular  to  the  xy-plane  which 
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contains  this  line  intersects  the  time  late  threat  density 
in  a  curve  whose  equation  is 


1       2 

1     c 


2    2^2   2   /  2~   2  2 

g^(y)  =  ( ^  e      ^       ^  )  (5^-J^^ ^ 

/t —  /   2    2   2  /2l\   m  a.,  a,, 

/27T  ^a   +  m  a^  X   Y 

2^22  2 

irr!l_!L!L!2L_w     ^  "y     .2. 

2^^2   2   2^^^"    2    2   2^^ 

m  a„  Qy         ay  +  ni  a„ 


2    2 
360      ^  (s-t,.)  cos  (e-tjj) -2s  .  t^-ycos  (9-i|;) 

I    exp[-  i   ^  ^  ^  ^ 


n  2  ^2   2 

0  m  a^ 


2    2 

(s.t-.)  sin  (6-'^) -2s  .  t^-ysin  (e-(|»)     ,, 

2  '  ^    360    ^    ' 

v—  C 

found  by  making  the  substitution  x  =  ^ m  the  density 

function  f    (x,y;t  )  and  expanding. 

X  /  X       "     Li 

For  computational  purposes,  assume  that  the  width  W, 
of  the  region  of  the  time  late  threat  ellipse  cut  out  by  the 
angular  subinterval  around  Bi^  will  be  nearly  constant  through 
the  ellipse.   Let  the  points  of  intersection  of  the  line  of 
bearing  Q,    with  the  time  late  threat  ellipse  be  (X,  wY,,), 
{X,^,Y,^).    Then  approximate  the  time  late  density  at  speed 
s .  over  this  region  by  the  absolute  value  of  the  product  of 
W,  and  the  area  under  the  curve  g,.  between  Y,  ,  and  Y,  ^.   The 
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area  under  curve  g   can  be  approximated  as  follows.   Sub- 
divide the  interval  (^  wY  2)  into  n,  segments  of  length  h. 
Evaluate  g   at  the  midpoint  of  each  segment,  y^.   Note  that 
the  first  term  in  parenthesis  in  (4)  is  a  constant  equal  to 

m  times  the  value  of  the  density  of  a  normal  random  variable 

2    2   2 
with  mean  zero  and  variance  a   +m  a„   evaluated  at  c.   The 

X  A 

second  term  in  parenthesis  is  the  density  function  of  a  normal 
random  variable  with  mean 

2 

ca^ 


and  variance 


2^22 
Oy   +  m  Oj^ 


2   2   2 

m  a^  a^ 

r~    2    2 

^Y   +  '^  <^x 


evaluated  at  y.   The  variable  y  also  appears  in  the  integral 

term  in  equation  (4).   Numerically  evaluate  this  term  of 

(4)  with  y  =  y,-   Let  g-r(yT.)  be  the  result  of  this  computation 

Then,  evaluating  the  area  under  the  curve  g^  between  Y,  , 

and  Y,  _  is  equivalent  to  the  calculation: 


1  1/  a,,  +  m  a.. 

r         m f C .  V      Y X    ./,  ^   ^   /„  X 

all   /  2   2  2  ^'/  2,  2  2    -^^^7; *'^J'  '^''^' 

y^Y         "^  ^X  l/'^Y    "^  ^X 


where  <))(•)  is  the  density  function  of  a  standard  normal  random 
variable,  and 
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^   2       I      2^2   2 

^J   -   (^J  -    2^  2   2^  ^    ma^o^ ^  ' 

Oy  +  m  a^        X  Y 


An  analagous  situation  arises  if  the  substitution  Y  =  mX  +  c 
is  made  for  y  in  the  time  late  density  f    (x,y;t  ). 

X  /  Y  Li 

8.   The  value  of  the  density  over  the  subinterval  con- 
taining 6,  is  then  weighted  by  the  instantaneous  probability 
that  the  bearing  error  is  B,^  -  3: 


^k2 


^k 
Y 


9.   Steps  C.7  and  C.8  are  then  repeated  for  each  sub- 
interval  of  (3  ,6rJ  and  the  results  of  each  calculation  are 

Li    U 

summed. 

10.  The  uncertainty  in  sensor  location  is  accounted  for 
by  repeating  steps  C.4  through  C.9  with  the  assumption  the 
sensor  is  located  at  each  of  the  three  points  in  Figure  5, 
multiplying  by  the  probability  the  sensor  lies  in  that 
interval  and  summing  each  result. 

11.  The  result  of  calculations  in  steps  C.l  through 
C.IO  is  the  likelihood  that  the  threat  lies  on  bearing  3 
given  the  threat  was  originally  located  in  the  p,  x  100% 
confidence  ellipse,  6   e  (Bt/S„)/  the  sensor  is  located  in 
the  Pt  X  100%  circle  and  the  threat  speed  is  s. .   The  con- 
dition that  the  speed  is  s.  is  removed  by  repeating  the 
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calculations  C.l  through  C.IO  for  each  of  the  speeds  s., 
i  =  1,...,M  multiplying  the  result  by  the  probability  that 
the  speed  equals  s.  and  summing  all  M  results.   The  final 
result  is  the  likelihood  that  the  threat  lies  on  6  given 
it  was  originally  located  in  the  p,  x   100%  ellipse, 
Qrj,  e    {Q    ,e>    ),    and  the  sensor  is  in  the  p-  x  100%  circle  esti- 
mated as 


M  3   , 

I    P[S=s.]  I      I  f  '(3)  . 
i=l         j=l     -' 


36 


IV.   CONCLUSIONS 

Possible  applications  of  and  extensions  to  the  algorithne 
developed  in  Chapter  III  are  discussed  in  this  chapter. 

As  indicated  in  Chapter  I,  the  objective  of  this  paper 
has  been  to  develop  a  means  of  assessing  the  likelihood  a 
threat  whose  position  has  previously  been  estimated  lies  at 
a  given  measured  bearing  from  a  local  sensor.   The  procedures 
were  developed  with  a  view  towards  permitting  the  user  to 
make  comparisons  among  lines  of  bearing  measured  by  different 
sensors  or  among  conflicting  bearing  information  generated 
by  one  sensor.   The  approach  chosen  has  been  to  estimate  the 
likelihood  that  a  threat  lies  on  bearing  6  given  that  the 
threat  is  located  in  an  ellipse  of  specified  confidence 
p,  X  100%,  that  6  is  measured  with  p_  x  100%  accuracy,  and  that 
the  sensor  position  is  measured  with  p^  x  100%  accuracy.   The 
algorithm  to  calculate  the  likelihood  in  the  cases  where  the 
probability  distribution  of  the  target  can  be  assumed  to  be 
bivariate  normal  at  the  time  of  the  bearing  measurement  has 
been  implemented  on  the  TI-59  calculator.   The  cases  in  which 
this  program  applies  are  the  following:   (1)  when  the  time 
elapsed  since  generation  of  the  threat  error  ellipse  is  small 
enough  to  justify  using  the  original  estimate  of  the  ellipse; 
(2)  when  threat  course  and  speed  are  known,  in  which  case 
the  ellipse  center  is  translated  from  the  original  position 
according  to  the  course,  speed  and  elapsed  time  information; 
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(3)  when  the  threat  can  be  assximed  to  be  moving  about  in  a 
random  manner  over  a  significant  region  in  such  a  way  that 
the  ellipse  center  remains  unchanged,  but  the  x  and  y  vari- 
ances have  increased.   When  none  of  these  cases  hold,  but 
the  course  is  assumed  uniformly  distributed  over  (0'',360°) 
and  speed  has  a  discrete  distribution  over  a  finite  interval, 
the  threat  density  at  the  time  late  t   is  not  a  bivariate 

Li 

normal  density.   The  algorithm  applicable  in  this  case  has 
not  been  implemented,  but  is  described  in  some  detail  in 
Chapter  III  and  Appendix  B. 

Once  the  appropriate  computation  has  been  completed  for 
each  of  the  bearings  considered,  the  results  can  be  used  to 
weight  the  value  of  several  bearings  in  refining  the  threat 
location  estimate  provided  by  the  external  sensor.   Note  that, 
although  unlikely,  f'(3)  may  correctly  be  greater  than  one. 
Comparisons  using  these  likelihoods  should  be  made  only  when 
the  upper  and  lower  bounds  on  the  true  bearing  fan  for  each 
bearing  are  chosen  at  the  same  probability  level  p_  and  the 
uncertainty  areas  for  all  sensors  include  the  same  probability 
level  p^.   Further,  these  likelihood  levels  should  be  selected 
so  as  not  to  exclude  a  significant  portion  of  the  appropriate 
density.   If  p^  or  p.  are  not  the  same  in  all  cases  to  be 
compared,  f'(6)  must  be  divided  by  the  applicable  value  of 
p   or  P-.  for  each  measured  bearing  3  to  be  considered. 

Having  established  the  relative  value  of  available 
bearing  information,  the  user  can  allocate  weapons  or  further 
search  effort  accordingly.   However,  the  probabilities  calculated 
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are  strictly  ordinal  data  and  do  not  define  a  redistribution 
of  target  location  probability  based  on  additional  informa- 
tion.  Further,  the  threat  ellipse  does  not  contain  the 
target  with  certainty.   The  power  to  predict  the  probability 
of  success  of  a  search  or  weapons  allocation  plan  based  on 
the  priorities  established  by  these  procedures  is  limited 
by  these  constraints.   In  this  area  in  particular  further 
research  would  be  useful. 

In  situations  such  as  that  for  which  the  procedures  in 
this  paper  have  been  developed,  where  a  track  has  not  been 
developed  on  the  target,  introduction  of  unrelated  bearing 
data  to  a  target  motion  model  could  impact  significantly  on 
the  reliability  of  future  position  predictions.   If  there  is 
high  confidence  in  the  reliability  of  the  estimate  of  the 
threat  ellipse  provided  by  the  ocean  surveillance  sensor, 
the  prioritization  established  herein  could  be  used  to  pro- 
cess bearing  data  prior  to  input  to  a  target  motion  model. 
Using  a  pre-established  threshold,  only  those  bearings  which 
coincide  with  the  threat  with  an  acceptable  level  of  like- 
lihood could  be  used  to  refine  or  update  a  track  on  the 
threat. 

Desirable  enhancements  to  the  algorithms  include  providing 
for  the  instances  in  which  the  interval  of  uncertainty  of 
the  target  course  is  known  to  be  less  than  (0°,360°).   Further, 
if  the  circular  region  of  radius  R  contains  the  sensor  with 
significantly  less  than  86%  confidence  or  the  assumption  of  a 
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bivariate  normal  distribution  of  sensor  location  is  unsatis- 
factory, it  is  left  to  the  user  to  modify  the  calculations 
accordingly. 

The  utility  of  these  algorithms  would  also  be  improved 
by  implementation  on  a  larger  and  faster  system  than  the 
TI-59  calculator. 
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APPENDIX  A.   TI-59  PROGRAM  VERBAL  FLOW  AND  USER'S  INSTRUCTIONS 


Part  I 

Step  Number  Verbal  Flow 


000  -  029  Enter  the  confidence  level  p,  for 

the  threat  ellipse.   Calculate  the 
value  of  k  for  the  given  p,: 


k  =   /-2  ln(l-p,)  . 


020  -  029  Enter  the  length  of  the  semi-major 

axis,  A.   Calculate  a„  =  A/k. 


0  30  -  052  Enter  the  length  of  the  semi-major 

axis,  B.   Calculate  a^  =  B/k. 

2        2 
Calculate  B/A,  o^a^,    o^     and  a„  . 


053  -  058  Enter  orientation  of  semi-major 

axis,  6. 

059-063  Enter  bearing  from  sensor  to  center 

of  threat  ellipse,  a.   Calculate 

9-a. 

064  -  073  Enter  distance  r  from  sensor 

position  to  center  of  threat 
ellipse.   Determine  rectangular 
coordinates  of  sensor  position 
(U,V)  from  polar  coordinates 
(-r,9-a) . 

074  -  089  Store  the  constants  360,  /2T  . 

Initialize  register  35  to  0. 
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090  -  096  Enter  number  of  standard  devia- 

tions desired  for  bearing  fan,  k  . 

09  7  -  104  Enter  standard  deviation  of  bearing 

error,  a..   Calculate  k^a„. 


105  -  143  Enter  the  angular  stepsize  desired, 

AS,  for  incrementally  stepping 
through  {Q^,Q    ) .      Calculate  jAB  . 
Calculate  the  largest  number  n  of 
increments  of  size  At3  contained  in 

^r,^ n    ~   T^B   degrees.   Initialize 

6  6    2 

counter  00  to  n+1.   Save  n+1  in 
register  20.   Determine 

6A6   =   kgOg  -  jA6   -  nA6  , 

the  residual  increment.   Calculate 

|6AS  . 

144  -  146  Initialize  counter  01  to  2. 

Calculations  will  be  made  at  the 
midpoint  of  each  interval  from 
6  +  0-A6  to  6  +  k-a^,  then  at  the 
midpoint  of  each  interval  from 
6  -  -jAB  to  6  -  k  a   ,  and  finally 
at  6.   Counter  01  indicates  whether 
calculations  are  complete  on  both 
sides  of  B. 


14  7  -  151  Enter  bearing  measured  by  sensor. 


152  -  154  Enter  standard  deviation  of  sensor 

position  a. 
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155  Enter  index  of  the  sensor  position 

to  be  used  for  this  run. 

156  -  157  Coordinates  of  the  sensor  position 

are  selected  in  accordance  with 
run  number  entered  above  in 
Subroutine  sin. 

158  -  238  Determine  whether  a  bearing  6' 

parallel  to  either  axis  is  included 
in  the  fan  (3j^,By).   If  (Sj^/Sy) 
includes  a  bearing  parallel  to  the 
y-axis,  use  program  1  for  Part  II. 
The  appropriate  program  number  is 
displayed  in  calculator  display 
register.   If  {Q-,Q^)    does  not 
contain  a  bearing  parallel  to 
either  axis,  use  program  1. 

239  -  259  Subroutine  P/R. 

260  -  340  Subroutine  sin. 
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Part  II,  Programs  1  and  2 

Step  Number  Verbal  Flow 


000  -  003  Initiate  6' . 

004  -  00  7  If  the  last  angular  increment  on 

this  side  of  6  has  been  considered, 
go  to  step  513.   Otherwise 
continue. 

00  8  -  012  Decrement  counter  for  angular 

increment.   If  counter  =  0,  go  to 
021.   Otherwise  continue. 

013  -  020  Remove  flag  to  indicate  this  is 

not  the  last  angular  increment. 
Recall  A6/  the  input  angular  step- 
size.   Go  to  029. 

021  -  028  Set  flag  to  indicate  this  is  the 

last  angular  increment  on  this 
side.   Add  one-half  the  input 
angular  stepsize  A6  and  one-half 
the  residual  stepsize  5A6: 

jAB  +  y6AB. 

029  -  030  Increment  6 '  by  the  appropriate 

stepsize. 

031  -  045  Convert  6'  to  an  angle  between  0° 

and  360°. 

046  -  051  Calculate  9-6'.   Print  6-6'. 
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052  -  064  If  Ie-6'l  =  90°  or  =  270°,  go  to  236 

065  -  074  If  |e-6'|  =  0°  or  =  180°,  go  to  075 

Otherwise  go  to  133. 

075  -  086  If  the  absolute  value  of  the  y 

coordinate  of  sensor  position  is 
greater  than  the  length  of  the 
semi-minor  axis,  Q,    go  to  004.   In 
this  case  6 '  does  not  intersect  the 
error  ellipse.   Otherwise  continue. 

087  -  088  Set  flag  2  to  indicate  that  the 

bearing  6 '  is  parallel  to  the 
X-axis. 

089  -  111  Calculate  the  coordinates  of  the 

points  of  intersection  of  S'  with 
the  threat  ellipse.   The  y-coordi- 
nates  are  equal  to  V,  the 
y-coordinate  of  sensor  position, 
x-coordinates  are  determined  in 
Subroutine  y  ,   The  points  of 
intersection  are  symmetric  about 
the  y-axis.   Thus,  X,  -,  =  -X,  2- 
Store  the  smaller  x  value  in 
register  27,  the  larger  in  register 
28. 

112  -  132  Save  locations  of  X^^^,  X^^^ '  ^j^i ' 

a„/  and  a„.   Go  to  291. 

133  -  138  Remove  flags  2  and  3  to  indicate 

that  3 '  is  not  parallel  to  either 
the  X-  or  y-  axes . 
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139  -  145  Calculate  the  slope  of  B ' : 

m  =   tan  (e-B  ' )  . 

146  -  154  Calculate  the  y-intercept  of  B ' : 

c  =  V  -  mU. 

155  -  176  .    If  c^  >  A^m^  +  B^,  go  to  004.   In 

this  case,  B'  does  not  intersect 
the  error  ellipse.   Otherwise 
continue. 


177  -  235  Calculate  X,  -,  and  X,  2/  the  x- 

coordinates  of  the  points  of  inter- 
section of  8 '  with  the  ellipse; 


\l  = 

-mc   + 

1  A^™^ , 

b2    - 

c2 

-mc   - 

A^^    +    b2 
a2 

\2      " 

1  A^™^ . 

b2    - 

2 

c 

A^m^   +  B^ 

Calculate  the  y-coordinates  of  the 
points  of  intersection  of  B'  with 
the  ellipse,  ^y^i'    ^^2* 

^kl   =  ^\l  ^  ^ 

\2      =   ^^k2  ^  ^  • 
Go  to  299. 

236  -  247  If  |u|,  the  absolute  value  of  the 

x-coordinate  of  sensor  position, 
is  greater  than  the  length  of  the 
semi-major  axis.  A,  go  to  004.   B' 
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does  not  intersect  the  error 
ellipse  in  this  case.   Otherwise 
continue. 

248  -  249  Set  flag  3  to  indicate  that  8'  is 

parallel  to  the  y-axis. 

2  50  -  2  72  Calculate  the  coordinates  of  the 

points  of  intersection  of  B '  with 

the  error  ellipse.   The  x- 

coordinates  are  equal  to  U,  the 

x-coordinate  of  sensor  position. 

The  y-coordinates  are  determined 

in  Subroutine  y  .   The  points  of 

intersection  are  symmetric  about 

the  X-axis.   Thus,  Y,  ,  =  -Y,  „. 

kl     k2 

Store  the  smaller  y  value  in 
register  29,  the  larger  in  30. 

273  -  29  8  Save  the  locations  of  Y, ^ /  Y,^, 

299  -  32  8  Determine  the  width  of  the  strip 

around  6 ' ,  W: 

If  flag  0  set, 

W  =   (d,+d2)  tan(j6AB), 
Otherwise 

W  =   (d-L+d2)  tan  ^AB)  . 

d,  and  d»  are  the  distances  from 
(U,V) ,  the  sensor  position,  to  the 
intersection  points  (^vi'^vi^  ^^^ 
^^k2'^k2^'  ^sspectively.   Distances 
are  calculated  in  Subroutine  log. 
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Multiply  W  by  A3  if  flag  0  not 
set.   Otherwise,  multiply  by  6Af 
Save  result  in  register  37. 

329  -  334  If  6'  is  parallel  to  either  the 

X-  or  the  y-axis,  go  to  468. 
Otherwise  continue. 


When  using  program  1; 

335  -  380  Express  the  y-coordinates  of  the 

intersection  points  as  standard 
normal  random  variables,  Z,  and  Z, 


Co 


7    -   ^i     2  _^  ^2    2 

Z.   -        ^Y   -"  ^   ^X 
[ml  o^a^ 


^/   2~7   2   2 
1/0^      "''  ^  ^X 

i  =  1,2. 

381  -  391  Sort  the  values  Z.  in  descending 

order.   Let  Z- '  be  the  larger 
value.   Z, '  is  the  smaller. 

392  -  401  Calculate  ^{Z^')    -  ^{Z^'),    the 

probability  that  a  standard  normal 
random  variate  lies  between  Z, ' 
and  Z^' . 
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402  -  424  Multiply  the  results  of  steps 

392-401  by: 


/2     2   2      I      2 
a    +  m  a        yo 


■)^  Iml  X  w. 


2   2 

X       y-Y      ^   "^    ^X 


where  ^    is  the  standard  normal 
density  function. 

The  result  of  the  calculations  in  steps  335-424  is 


^k2 

^ki 


■X,Y   m 


When  using  program  2; 

335  -  372  Express  the  x-coordinates  of  the 

intersection  points  as  standard 
normal  random  variables  Z-,  and  Z^ 

2 

mca^ 

\2^ 


J      2~      2      2 

1/a    +  m  a 

Z.   = 

1 


^X'^Y 


,/ 2  ,   2   2 

U  a         +  m  a 


i  =  1,2 


373  -  383  Sort  the  values  of  Z.  in  descending 

order.   Let  Z'    =   larger  value. 
Z, '  =  smaller. 
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384  -  393 


Calculate  ^(Z^')  -  $(Z^'),  the 
probability  that  a  standard  normal 
random  variate  lies  between  Z,  ' 
and  Z'. 


394   -    413 


Multiply   the   results   of   steps 
384-393   by: 


*(■ 


J      2~~T2      2  ./      2~      2      2 


-)     X    W, 


where  (p    is  the  standard  normal 
density  function. 

The  result  of  the  calculations  in  steps  335-413  is 

Y, 
W 


•k2 


f     f^   (x,mx+c)  dx 

r  A  ,  X 

kl 


Regardless  of  which  program  is  in  use: 


425  -  465 


Multiply  results  of  previous 
calculation  by 


where  B " 


S' 


,    B'  <  180« 


B'-360°   ,  B'  >  180° 


and  (J)  is  the  standard  normal  density, 


The  result  of  this  calculation  is 


fg(B;B')   |w  /  fx^y(^'y)  ^y 
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466  -  467  Go  to  508 

46  8  -  50  7  Calculate  the  area  under  the  curve 

formed  by  the  intersection  of  6 ' 
with  f^  ^{x,y) .      When  3'  parallels 
the  X-axis  this  calculation 
becomes : 

^   4'(;7-)  [HZ  )  -  $(Z,)]  , 

Y       Y 

where  Z^   =  X^^/a^   and  Z^  =  X^^/a^    , 

V  is  the  y-coordinate  of  the  sensor 

position. 

When  3'  is  parallel  to  the  y-axis , 

the  calculation  is: 

^  ^{^)  [HZ^)    -    $(Z,)  ]  , 

where  Z^  =  Yj^2/'^y  ^^^  ^1  ^  ^kl'^'^Y  ' 

U  is  the  x-coordinate  of  the  sensor 

position. 

Go  to  422  to  complete  calculation 

of   fg(3;3')  jW  /f^^Y(x,y) 1  . 

(Go  to  410  in  Program  2) . 

50  8-512  Accumulate  the  probability  at  each 

angular  interval.   Display  result. 

513  -  515  If  flag  1  is  set  go  to  558. 

Otherwise  continue. 

516  -  520  If  flag  0  is  not  set,  that  is  if 

the  angular  increment  just  con- 
sidered was  not  the  last  on  this 
side  of  3,  go  to  004  and  continue 
calculation  on  same  side  of  3. 
Otherwise  continue. 
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521  -  524  Decrement  the  counter  01.   If 

counter  is  now  0  go  to  54  3.   In 
this  case,  the  probabilities  have 
been  calculated  on  both  sides  of  B. 
The  calculation  at  6  remains  to  be 
done.   Otherwise  continue. 

525  -  542  Remove  flag  0.   Multiply  A3,  ^A6, 

1  2 

and  jiSAB  by  -1.   Reinitialize 

counter  0  0  to  n+1.   Go  to  000  to 

begin  calculation  on  second  side 

of  S. 

543  -  567  Set  flag  01  indicating  that  calcu- 

lations on  both  sides  of  6  have 
been  completed.   The  next  iteration 
will  do  the  calculation  at  6'  =  B. 
Remove  flag  00.   Initialize  B'  to  0 
Recall  6.   Go  to  029. 

558  -  566  Remove  flag  01.   Display  the 

accumulated  likelihood.    STOP 


The  result  of  this  calculation  is  approximately  f.'(B)  with 
sensor  at  (U.,V.)  where  j  =  run  number  entered  in  Part  I. 


567  -  587  Subroutine  y^  calculates  the  points 

of  intersection  when  Q '    is  parallel 
to  either  the  x  or  y  axis: 


Xki  =  -^  z/i  -  <& 


\2  =  *  i/i  -  <!'' 
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Parallel    to   y-axis: 


\l    =    \2    =    " 


-.1  =  -B  iT"^ 


Y,.,  =  s/rr^ 


kl 


589  -  610  Subroutine  log  calculates  the 

distance  from  sensor  position  to 
point  of  intersection  of  6 '  with 
the  ellipse: 


I.   =  t/Tu-X,  .  )  ^  +  (V-Y,  .) 
1     ^  ki         ki 
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part  III 

Step  Number  Verbal  Flow 


000  -  012  Enter  result  of  run  1. 

Multiply  by  1/3. 

013  -  021  Enter  result  of  run  2. 

Multiply  by  1/3. 

022  -  029  Enter  result  of  run  3. 

Multiply  by  1/3. 

030  -  032  Display  likelihood. 

033  -  041  If  p_  is  the  same  for  all  bearings 

to  be  compared,  enter  1.   Go  to 
042.   Otherwise,  enter  p_  for  this 
bearing.   Divide  likelihood  by  p^. 
Display  result. 

042  -  050  If  p^  is  the  same  for  all  bearings 

to  be  compared,  enter  1  and  STOP. 
Otherwise,  enter  p^  for  this  bearing. 
Divide  likelihood  by  p^.   Display 
result.   STOP. 
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USER'S  INSTRUCTIONS 

The  program  to  determine  the  likelihood  that  the  threat 
lies  along  bearing  6  given  the  threat  is  in  the  confidence 
ellipse.  Sly,  £  (6^,3^)  and  the  sensor  is  within  the  p  x  100% 
confidence  circle  is  in  three  parts.   All  parts  require  the 
use  of  a  printer  and  the  use  of  the  Applied  Statistics 
Library  Module.   Prior  to  running,  the  calculator  must  be 
repartitioned : 

1.  Enter  4 

2.  Press  2nd  OP  17 

Part  I 

1.  Read  sides  1  and  2  of  Part  I 

2.  Read  side  4  of  Part  II  either  program  1  or  program 

2.   Since  program  2  is  used  more  often,  unless  it  is  known 

that  the  bearing  fan  (3,6)  contains  a  bearing  parallel  to 

J-i   u 

the  major  axis,  recommend  using  program  2  of  Part  II. 

3.  Press  RST 

4.  Enter  p, ,  the  confidence  level  of  the  threat  ellipse, 
Press  A 

5.  Enter  A,  length  of  the  semi -major  axis.  If  data 
provided  is  length  of  the  entire  major  axis,  divide  by  2 
before  entering. 

Press  B. 

6.  Enter  B,  length  of  the  semi-minor  axis.   If  data 
provided  is  length  of  entire  minor  axis,  divide  by  2  before 
entering.   Press  C. 
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7.  Enter  e,  the  bearing  of  major  axis  from  North. 
0°  <_  e  _<  180".   e  is  entered  in  degrees.   Press  D. 

8.  Enter  a,  the  bearing  of  the  threat  ellipse  center 
from  the  estimated  sensor  position.  ^°   ±  ol    <_   360",  in 
degrees.   Press  E. 

9.  Enter  r,  the  range  from  the  estimated  sensor  posi- 
tion to  the  center  of  the  threat  ellipse.   Press  2nd  A'. 

10.  Enter  k    ,    the  number  of  standard  deviations  desired 

p 

in  one  direction  from  s.   The  program  will  construct  a  fan 
of  equal  size  on  the  other  side  of  6-   Press  2nd  B'. 

11.  Enter  a  ,  the  standard  deviation  of  bearing  error. 
Press  2nd  C . 

12.  Enter  ^s,,    the  desired  angular  stepsize.   Press  2nd 
D'  . 

13.  Enter  Q,    the  measured  bearing.   Press  2nd  E'. 

14.  Enter  a,    the  standard  deviation  of  sensor  position. 
Press  R/S. 

15.  Enter  the  number  of  this  run,  1,  2  or  3.   When  run 

number  =  1,  sensor  is  located  at  its  estimated  position 

(u_,v-,).   When  run  number  =  2,    the  location  will  be 
o   o 

(u+ .97acos (9-8-90) ,v  + .97asin (9-6-90) ) .   When  run  number  =  3, 
the  location  will  be  (Uq- . 97acos (9-6-90) ,Vo-.  97asin (9-6-90) ) . 
Press  R/S.   If  program  number  displayed  matches  that  of  the 
Part  II  side  4  read  in,  continue  to  16.   Otherwise,  press 
2nd  CMS,  RST.   Read  in  side  4  of  the  Part  II  program  which 
corresponds  to  the  number  displayed.   Repeat  3  through  15. 
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16.   Press  4  2nd  WRITE.   Rerecord  side  4  of  the  Part  II 
program  read  in.   This  enables  data  entered  in  Part  I  to 
be  transferred  to  Part  II. 

Part  II 

1.  Read  sides  1,  2,  3  and  4  of  Part  II  Program  1  or  2 
as  selected  by  the  Part  I  program.. 

2.  Press  RST 

3.  Press  R/S.   Values  of  B.  and  w|f„(S;6^)f^     (6^) 

^       ^  6     t   B^j^uE   ^ 
will  print  alternately.   Final  result  also  prints  out  at  end 

of  calculation. 

4.  Record  final  result:   Likelihood  threat  lies  along 
bearing  6  given  threat  is  in  p,  x  100%  ellipse,  true  bearing 
is  in  Pp  X  100%  fan  and  sensor  is  at  location  used  in  this  run, 
Parts  I  and  II  must  be  completed  3  times  (Run  numbers  1,  2 
and  3)  before  proceeding  to  Part  III  if  uncertainty  in  sensor 
position  is  being  considered. 

Part  III 

1.  Read  side  1  of  Part  III. 

2.  Enter  result  of  run  1.   Press  A. 

3.  Enter  result  of  run  2.   Press  B. 

4.  Enter  result  of  run  3.   Press  C. 

5.  Enter  p^ ,  confidence  level  of  {&,$),    if  necessary. 

Z  Li  U 

Otherwise  enter  1.   Press  D. 

6.  Enter  p.  confidence  level  of  sensor  position,  if 
necessary.   Otherwise,  enter  1.   Press  E. 
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APPENDIX  B.   TIME  LATE  VERBAL  FLOW 


VERBAL  FLOW 


1.  Enter  p,  .   Determine  the  value  of  k:   k  =  i/-2  In  (1  -  p,  ) 

2        2 

2.  Enter  A.   Determine  a^  =  (A/k) 

2        2 

3.  Enter  B.   Determine  a^   =  (B/k) 

4.  Enter  9,  orientation  of  major  axis 

5.  Enter  a,  bearing  from  sensor  to  ellipse  center 

6.  Enter  range  r  from  sensor  to  ellipse  center 

7.  Calculate  (u  ,v  )  the  coordinates  of  mean  sensor 

o   o 

position.   u  =  -r  cos(9-a),   v  =  -r  sin(0-a) 

8.  Enter  a^,    the  standard  deviation  of  bearing  error 

p 

9.  Enter  k  ,  the  number  of  standard  deviations  to  be 
included  in  the  bearing  fan  on  each  side  of  g. 

10.  Enter  desired  angular  stepsize,  A6 

11.  Calculate  the  number  of  iterations  of  size  A6  required 
on  each  side  of  s :   I  =  [  {k   a  .  -  ^/\Q)  /  AQ]    where  [•]  means 

the  greatest  integer  less  than  or  equal  to  the  value  within. 

12.  In  general  (k  a^  -•5-AB)/A3  is  not  an  integer.   Determine 

p  p   ^ 

the  size  of  the  fractional  increment: 

5A3  =   ( (k^a  -  2A6)/A6  -  I) A6. 

13.  Enter  6,  the  measured  bearing.   Let  PROB2=0,  PROB1=0 

14.  Enter  the  time  late  t^. 

15.  For  each  of  the  discrete  threat  speeds  to  be  con- 
sidered, repeat  steps  16  to  43.   Then  go  to  44. 

16.  Enter  the  target  speed  s.  Let  PROB  =  0 
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17.  Expand  the  outer  limit  of  the  threat  ellipse: 

Let   A'   =   A  +  st^ 
B'   =   B  +  St^ 

Li 

18.  Let  6'  =  3 

19.  For  each  of  the  I  increments  of  size  as  repeat 
steps  20  to  38. 

20.  Let  6'  =  6'  +  A6 

21.  Calculate  (e-S').   If  |6-6'|  =  0°  or  180°,  go  to  27, 
If  |e-e'|  =  90°  or  270°,  go  to  30 

22.  Calculate  m  =  tan  (e-6'),  the  slope  of  the  line  of 
bearing. 

23.  Calculate  c  =  mU+V,  where  (U,V)  is  the  sensor  loca- 
tion for  this  iteration.   c  is  the  y-intercept  of  the  line 

of  bearing. 

2  2    2    2 

24.  Calculate  m  A  +B  -  c  .   If  this  quantity  is  less 

than  zero,  the  line  of  bearing  s'  does  not  intersect  the 
threat  ellipse.   Go  to  step  20  and  process  next  increment  of 
size  A3,  if  any  remain.   If  all  I  intervals  of  size  A3  on 
this  side  of  3  have  been  considered  go  to  step  40. 

25.  Calculate  the  points  of  intersection  of  the  line 
of  bearing  with  the  expanded  ellipse: 

(-mc  +  (B'/A')/A'^m^  +  B'^  -  c^)  (A'^) 
1  A'  m   +  B ' 


(-mc  -  (B'/A')/A'^m^  +  B ' ^  -  c  )  (A'^) 
^2   =  A- V  +  B'^ 
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Y^   =   mX,  +  c,    Y2   =   mX^  +  c 

26.  Go  to  32 

27.  If  |v|  >  B",  the  bearing  6'  parallel  to  the  x-axis 
does  not  intersect  the  threat  ellipse.   Go  to  step  20  and 
process  next  increment  of  size  A6  if  any  remain.   If  all  I 
intervals  if  size  A3  on  this  side  of  3  have  been  considered, 
go  to  40. 

28.  Calculate  the  points  of  intersection 


Y,   =   Y2   =   V 


:^   =  -B'jjl  -    (V^/B'^) 


X2   =   B'/l  -  (V^/B'^) 

29.  Go  to  32 

30.  If  |ul  >  A',  the  bearing  3'  parallel  to  the  y-axis 
does  not  intersect  the  threat  ellipse.   Go  to  step  20  and 
process  the  next  increment  of  size  A3  if  any  remain.   If 
all  I  intervals  of  size  A3  on  this  side  of  3  have  been  con- 
sidered, go  to  40. 

31.  Calculate  the  points  of  intersection: 


X^   =   X2   =   U 


'^       =   -A'  t/l  -  (U^/A'^) 


Y^  =  A'  /i-  (uVa'^; 
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32.   Calculate  the  median  width  of  the  strip  of  the 
ellipse  defined  by  the  angular  subinterval  under  considera- 
tion.  If  the  subinterval  is  of  size  A6: 

d  +d 
W  =   -="2 — -{2    tan  |-A6) 

If  the  subinterval  is  of  size  5A3: 

d^  +  d2       1 
W  =  -^^ — -{2    tan  j6A6) 

In  these  expressions  d-  is  the  distance  from  the  sensor  to 


the  ith  point  of  intersection,  i  =  1,  2:   d.  =]/(U-X.)  +  (V-Y . ) ^ 

33.  If  |e-3'|  =  0"  or  180°,  go  to  36.   If  |  9-e  '  I  =  90° 
or  270°,  go  to  35. 

34.  Evaluate  the  target  density  from  Y, '  to  Y^ '  along 
the  line  y  =  mx  +  c: 

Let  Y^'   =  min(Y^,Y2) 

Y^'   =   max(Y^,Y2) 

(l)(z)   =   density  function  of  a  standard  normal 
random  variable  evaluated  at  z 

/'  2    2   2       C     2    2   2 
a^  +  m  a^  )  ^[c/ijo^     +  m  a^  ) 

Subdivide  the  interval  (Y  ',Y  ')  into  n  segments  of  length 
h.   At  the  midpoint  of  each  segement,  y.,  compute: 
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^     2 

r~2       2     2  Y  /     2       2     2 

hK(7/aY    +m  Gj^    71110^0^)    ^((y^ 2 — ~2 — 2^  ^^°Y    "''"'■  '^X   ^'^'^^^^ 


360            ,     (st. )    CX3S    (e-i^j)  -  2s1:L.y .  cxDs(9-ij;) 
/    e^[-7( 2—2 ^ 


m  o^ 


2       2 

(St.)    sin   (6-ij;)  -  2s-ty.   sin(0-ij;)        ,, 

2  ^^360 


°Y 


when  the  integral  must  be  numerically  evaluated. 
Sum  the  results  over  all  n  segments. 
Go  to  37. 

35.   Evaluate  the  target  density  from  Y, '  to  Y^ '  along 
the  line  x  =  U: 

Let  Y-^'  =  min(Y^,Y  ) 
Y2'  =  max(Y^,Y2) 
(|)  (z)  =  density  function  of  standard  normal  random 

variable 
K   =   -L  ,(iL) 

°x   ^x 

Subdivide  the  interval  (Y, ',Y2')  into  n  segments  of  length 
h.   At  the  midpoint  of  each  segment,  y.,  compute: 

T   y.    ^^°     T  (st-)^cos(e-ij;) -2st^Ucos(e-iij) 
hK— (*(-J-))   ]■   exp[-^( ^ ^ ^ 


^X   '^Y    .  a. 


av   a, 

0  "X 

2   2 

(st^)  sin  (e-ij;)  -  2st^y.sin(9-ij^)    , 

"^  '  2  ^^  360 
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where  the  integral  must  be  numerically  evaluated. 
Sum  the  results  over  all  n  segments.   Go  to  37. 

36.   Evaluate  the  target  density  from  X,'  to  X- '  along 
the  line  y  =  V: 

Let  X^'  =  min(X^,X2) 
X2'  =  max(X^,X  ) 
(j)(z)  =  density  of  standard  normal  random 
variable 
K   =   i  *(X, 


Subdivide  the  interval  (X, ',X  ')  into  n  segments  of  length 
h.   At  the  midpoint  of  each  segment,  x . ,  compute: 

2 

,    X.     360      ,  (st.)  cos(e-t)  -  2st.x.  (cos(e-ii;) 

hK—  (^(-J-))   f    exp[-^( 5 ^^ 

^X    ^X    6         2  aJ^ 

(stL)^sin^(9-'J;)  -2st^Vsin(9-iJ;)    ^ 
+  2  ^  ^  360 


where  the  integral  must  be  numerically  evaluated.   Sum 
results  over  all  n  segments. 

37.  Multiply  the  value  of  the  target  density  just 
computed  by  W. 

38.  Multiply  this  result  by 


^    (■ 


o  „  a , 


where  <P    is  the  standard  normal  density  function. 
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39.  Let  PROB  =  PROB+  (the  results  of  the  calculations 
in  steps  20  through  38  for  each  if  the  I  subintervals  of 
size  A6) . 

40.  If  the  fractional  interval  of  size  6A6  on  this 
side  has  been  considered  go  to  41.   Otherwise,  let 

B'  =  j6AQ    +   -jAS.   Repeat  calculations  21  through  38  once. 
Let  PROB  =  PROB+  (the  result  calculated  at  this  step) . 

41.  If  the  computations  on  both  sides  of  6  have  been 
computed,  go  to  42.   Otherwise  repeat  the  computations 
from  18  to  40  on  the  other  side  of  B  by  letting  ab  =  -AB/ 
5AB  =  -6AB. 

42.  Let  Q'    =    Q.      Repeat  steps  21  through  38  once. 
Let  PROB  =  PROB+  (the  results  of  this  calculation) . 

43.  Let  j  =  the  number  of  discrete  threat  speeds  to 
be  considered.   Let  PROBl  =  PROBl  +  P[S=s]PROB.   Go  to  15. 

44.  Repeat  steps  15  through  43  once  with  the  sensor 
located  at  each  of  three  points: 

(1)  (U,V)  =  (u  ,v  ) ,  the  mean  of  the  sensor 

o   o 

density; 

(2)  (U,V)  =  (u^ +  .97aCOS(9-B-90) , 

V  +  .97asin(e-B-90) ) ; 

(3)  (U,V)  =  (u  -.97  cos(e-B-90), 

V  -  .97  sin(9-B-90) ) ; 
o 

where  a  is  the  standard  deviation  of  the  sensor  density. 
The  value  of  PROBl  calculated  at  each  iteration  will  be 
weighted  by  the  approximate  probability  that  the  sensor  is 
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located  in  the  region  of  the  sensor  error  circle  repre- 
sented by  the  applicable  value  (U,V) .   If  p_  is  of  the 
order  .86  or  greater,  multiply  by  1/3. 

45.   The  value  of  PR0B2  calculated  after  completion 
of  step  44  is  the  relative  probability  that  the  threat 
lies  along  bearing  6  at  time  t  ,    given  the  threat  was  in 

Li 

p,  X  100%  ellipse  at  time  t  ,  6.  is  in  (3  ,6  )  and  the 
sensor  is  in  a  p^  x  100%  confidence  circle. 
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